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Summary. The construction of effective Hamiltonians describing corrections to 
fiat space particle dynamics arising from the granularity of space at very short 
distances is discussed in the framework of an heuristic approach to the semiclassi- 
cal limit of loop quantum gravity. After some general motivation of the subject, a 
brief non-specialist introduction to the basic tools employed in the loop approach 
is presented. The heuristical semiclassical limit is subsequently defined and the ap- 
plication to the case of photons and spin 1/2 fermions is described. The resulting 
modified Maxwell and Dirac Hamiltonians, leading in particular to Planck scale 
corrections in the energy-momentum relations, are presented. Alternative inter- 
pretations of the results and their limitations, together with other approaches are 
briefly discussed along the text. Three topics related to the above methods are 
reviewed: (1) The determination of bounds to the Lorentz violating parameters in 
the fermionic sector, obtained from clock comparison experiments. (2) The calcu- 
lation of radiative corrections in preferred frames associated to space granularity 
in the framework of a Yukawa model for the interactions and (3) The calculation 
of synchrotron radiation in the framework of the Myers- Pospelov effective theories 
describing Lorentz invariance violations, as well as a generalized approach to ra- 
diation in Planck scale modified electrodynamics. The above exploratory results 
show that quantum gravity phenomenology provides observational guidance in the 
construction of quantum gravity theories and opens up the possibility of probing 
Planck scale physics. 



1 Introduction 

Most theories of gravity suggest that our notion of space-time as a contin- 
uum needs to be revised at short distances (high energies) of the order of the 
Planck length l P sa 1CT 33 cm (Planck mass, M P « 10 19 GeV). At these scales 
quantum effects should be important and since space-time is to be consid- 
ered as a set of dynamical interacting variables, instead of a mere background 
where physics occurs, the quantum nature of them together with their cor- 
responding fluctuations could induce modifications to our standard notion of 
space-time. The consequences of space-time been considered as a continuum 
have been already successfully probed up to the much lower energies of « 10 3 
GeV, corresponding to the standard model of particle physics. The sixteen 
orders of magnitude between our current experimental arena and the region 
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where such new effects woufd become relevant provides in fact plenty of room 
to look for modifications of our generally accepted ideas of space-time. We 
will refer to this new possibility as space-time presenting a granular or foamy 
structure in the following. Central to this question is the long time honored 
problem of finding a consistent unification of gravity and quantum mechan- 
ics, which nowadays is been actively pursued by loop quantum gravity 
and string theory among others lines of research. At any rate, a com- 
plete formulation of quantum gravity must explain how the standard notion 
of space-time at macroscopical length scales is recovered, thus validating the 
many successful test of classical Einstein gravity that have been performed 
and which have served to make sure that the correct starting point to con- 
struct the full theory has been taken. 

Then a reasonable question to ask is whether or not such short length 
(high energy) effects leave any imprint in the dynamics of particles at stan- 
dard model energies, which we could be able to detect with present day 
technology and observational sensibilities. 

On one hand, from a purely phenomenological point of view one would 
expect some modifications to arise in the same way as particle propagation 
properties change when they move in a medium, with respect to those in the 
vacuum. Of course the analogy of a modified structure of space-time with a 
propagating medium is at most very tentative because the idea of a medium 
presupposes the existence of something external in which it is embedded. On 
the contrary, space-time is the arena where phenomena occur and there is 
nothing external to it. Nevertheless it is very plausible that a drastic change 
in our description of space-time would induce modifications in the way we 
deal with physics at such very short scales. Here we take the point of view 
that there is a remnant of such modifications at standard model energies, 
which are described by an effective field theory valid up to scales much lower 
than the Planck mass. These corrections manifest themselves as additional 
terms contributing to the propagation and interactions of known particles. 

On the other hand, for a long time it has been taken for granted that 
Planck scale phenomena is completely out of reach from present experiments 
and/or observations (astrophysical, for example). Recent investigations show 
that this is definitely not the case, thus opening the door to a new area of 
research called Quantum Gravity Phenomenology (QGP) that is designed to 
use existing and forthcoming experiments and/or observations to restrict or 
constraint competing theories of quantum gravity based on their predicted 
imprints at these lower energies [3]. A partial list of references describing 
such efforts is given in Ref . 0j |HJ ■ This is definitely a great advance over 
the purely aesthetically criteria that prevailed before. Even though QGP has 
made independent advances in restricting the different parameters encoding 
some of the proposed modifications induced at standard model energies, there 
still remain the open problem of calculating them as a rigorous semiclassical 
approximation of a fundamental theory. Only after this gap is filled one could 
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really be in position to use the observations to restrict the latter. The work 
presented here summarizes an heuristical step towards the estimation of the 
flat space semiclassical limit in the realm of LQG, together with further 
elaborations upon the results obtained in this way. Also some closely related 
topics originating from different approaches are discussed. 

Before going to the details, let us briefly describe some of the alternative 
points of view to be found in the literature regarding the question of whether 
or not a granular or foamy structure of space would induce modifications to 
particle dynamics at standard model energies. Such modifications have been 
mainly understood in terms of the violation of standard Lorentz covariance, 
either through the introduction of preferred reference frames or by means of 
the inclusion of an extended or deformed Lorentz relativity. In the following 
we will generically refer to theses possibility as Lorentz invariance violation 
(LIV). 

A first viewpoint is that modifications that manifest themselves in the 
form of LIV, do not arise. In this way covariance under the standard Lorentz 
transformations would be perfectly compatible with a discrete nature of space 
at Planck scale, which is encoded in the discrete spectrum of area and volume 
operators in LGQ, for example. This would be analogous to the well known 
property that the discrete spectrum of the angular momentum operator does 
not preclude the invariance of a system under the continuous rotation group 

A second possibility is that corrections in the form of LIV do in fact arise. 
Within this point of view different alternatives have been also considered, 
which can be separated in two further categories: 

(i) The first one can be characterized by a phenomenological parameter- 
ization of all possible corrections terms, according to the dimensionality of 
the corresponding LIV operators, which are assumed to arise via a sponta- 
neous Lorentz symmetry breaking of a more fundamental model like string 
theory, for example These vacuum expectation values define a set of 
preferred frames, called concordant frames UJ, in which the LIV terms can 
be maintained appropriately small when going from one frame to another 
via a passive (observer) Lorentz transformation. This is analogous to the 
description of atomic phenomena in the presence of an external magnetic 
field, where rotational invariance is broken by active (particle) transforma- 
tion. Nevertheless, one can rotate the apparatus and perform the experiment 
in the presence of the rotated external field. The Standard Model Extension 
(SME) JI] belongs to this class and it has been highly successful in allowing 
a unified description of the great amount of experimental data that has been 
gathered since 1960 ^H] m relation to the experimental verification of the 
isotropy of apace, the transformation properties among inertial frames and 
the validity of the discrete transformations C, P and T as well as its possi- 
ble combinations • The SME has been recently generalized to incorporate 
gravity ^Sj- More recently LIV models based upon dimension five operators 



4 Luis F. Urrutia 



have been constructed and thoroughly analyzed in Refs. ^1 ■ Within 
the approach (i) we find also direct extensions of Dirac and Maxwell equa- 
tions incorporating modifications which go beyond effective field theories ^7] . 
An alternative approach to LIV via spontaneous symmetry breaking can be 
found in Ref. [IB]. 

(ii) A second broad category of researches which leads to dynamical cor- 
rections can be distinguished by an effort to obtain such modifications either 
from an extension or deformation of the standard Lorentz relativity principle 
|19l 121)1 |2"T1 1221 |2"51 I24| . or from a fundamental description of quantum grav- 
ity like, for example: an effective field theory description of quantum general 
relativity [213 EH], a version of non- critical string theory [231211 EH or the 
approach of loop quantum gravity [123 EH1 021 CHU 133 EH] ■ Alternative ways 
of incorporating such dynamical modifications can be found in Ref. |3f>l \'M\ . 

The most immediate way in which the corrections to the dynamics show 
up is through modifications to the energy-momentum relations for the parti- 
cles EHUm 
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Eqg 

as shown above for photons and fermions, respectively. The connection of 
such modifications with gravity, together with the possible astrophysical ob- 
servation of the energy-dependent velocity in high-energy photons arriving 
from cosmological sources was suggested in Rcf. 40 . In fact, the frequency 
dependent photon velocity 

( 3 Ifel \ 

predicts a time delay At in the arrival of two photons having an energy 
difference AE and travelling a distance L, given by 

Eqg c 

Here Eqg denotes the scale where quantum gravity effects become rele- 
vant, which is usually taken as the Planck energy. As emphasized in Ref. 
|4"0"] . the tiny contribution of £,/Eqg can be amplified by selecting large en- 
ergy differences or, more effectively, distances L of cosmological magnitude. 
It is important to emphasize that this time delay is a purely kinematical 
consequence of the modified dispersion relations. The orders of magnitude 
L w W w l.y., AE « 20MeV, E QG = 10 19 Gey, £ « 1, lead to At « 10" 3 s, 
which is within the range of sensitivities St in actual and forthcoming gamma 
ray bursts (GRB) observations. In order to measure such effect it is necessary 
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Table 1. Bounds upon the quantum gravity scale Eqg/£ 



z 


6t(s) 


Objects 


E QG /i{GeV) 


Source 


0.031 


280 


Markarian421 


> 4.0 x 10 16 


Ref.|U 


0.0085-3.42 


0.064 


GRB's (BATSE) 


> 6.9 x 10 15 


Ref.gT] 


0.3 


0.031 


GRB021206 (RHESSI) 


> 1.3 x 10 17 


Ref.gg 



that St < At. Present sensitivities allow for estimates of the lower bounds for 
Eqg/£, presented in Tabled Future planned observations of GRB at cosmo- 
logical distances having well determined red shifts z, together with greater 
sensitivities ranging form 1CP 6 s (RHESSI: Reuben Ramaty High Energy So- 
lar Spectrometer) to 10~ 7 s (GLAST: Gamma Ray Large Area Telescope ) 
will allow a substantial increase of such bounds. The fireball model of GRB 
emission 021 HH predicts also the generation of 10 5 — 10 10 GeV neutrino 
bursts which will be detected by observatories like NUBE (Neutrino Burst 
Experiment), OWL (Orbiting Wide-angle Light collector experiment) and 
EUSO (Extreme Universe Space Observatory), for example. This will open 
up the possibility of using such particles, some times detected in coincidence 
with the respective photons, to set further bounds upon the quantum gravity 
scale Eqg- 

Very soon after the proposal of Ref.0O]> the first derivation of a consis- 
tent electrodynamics leading to the dispersion relations was obtained in 
the framework of a loop quantum gravity inspired model |3U|. The resulting 
modified Maxwell equations are 

V-£7 = 0, d t E = -V x B + 2££ P \7 2 B, (5) 
V-B = 0, d t B = +V x E -2££ P \7 2 E, (6) 

leading to the energy-momentum relation 

u,± = |fc|(lT2^p|fc|). (7) 

Additional bounds upon the parameters describing the quantum grav- 
ity induced modifications have been obtained by incorporating the dynamics 
through the SME or similar constructions. In particular, the topics of Lorentz 
and CPT violations have been thoroughly studied in low energy physics via 
theory and experiments related to: Penning traps, clock comparison measure- 
ments, hydrogen- antihydrogen studies, spin polarized dispersion and muon 
experiments, among others subjects |4~K] . 

Finally we mention the use of astrophysical phenomena to discuss such 
modified theories. Distinguished examples are the bounds imposed by po- 
larization measurements from astrophysical sources |46| . the study of ultra 
high energy physics processes, among them cosmic rays, 03 0BI B2I 0231 E] 
and the consequences of the detected synchrotron radiation from the Crab 
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nebulae, as well as that from other objects |15U16ll5*21l53) . For recent reviews 
about such topics see Ref.fS"!]. 

This contribution summarizes the work carried over in collaborations with 
different colleagues and it is organized as follows. Section [3 contains a non- 
specialist review of some basic elements of LQG to be subsequently employed 
in the estimation of the induced dynamical modifications. In section we de- 
scribe the general features involved in the heuristical calculation of the effec- 
tive photon and fermion Hamiltonians incorporating Planck scale corrections 
induced by LQG. These Hamiltonians are presented in subsection 3.3 to- 
gether with some related comments. The remaining Section 0] is devoted to 
selected phenomenological applications of models presenting Planck scale dy- 
namical corrections. Using existing data from clock-comparison experiments, 
in subsection 4.1 we obtain stringent bounds upon combinations of param- 
eters appearing in the effective fermionic Hamiltonian previously derived. 
The incorporation of radiative corrections to the description of LIV effects 
in preferred frames associated to space granularity is discussed in subsection 
4.2, leading to severe fine-tuning and naturalness problems. Subsection 4.3 
contains the discussion of synchrotron radiation in LIV electrodynamics, em- 
phasizing the model of Ref. ^3] , which is phenomenologically formulated as 
a theory parameterizing LIV with dimension five operators. Finally, a unified 
description of radiation in Planck scale modified electrodynamics including 
different models is presented in subsection 14.41 

2 Basic elements from Loop Quantum Gravity (LQG) 

Detailed reviews of LQG can be found in Refs.(I]|5H|. Here we just indicate 
some of the basic features that will be relevant in our heuristic calculation and 
present an intuitive, non-specialist introduction to this approach. LQG turns 
out to be the formulation of Einstein gravity as a Hamiltonian gauge theory 
with additional constraints, written in terms of non-local gauge covariant, 
diffeomorphisms invariant quantities. It is formulated in four dimensions and 
the matter couplings are obtained by rewriting the standard ones in terms 
of the new variables. This background independent, non-perturbative theory 
has allowed substantial progress in the old problem of producing a consistent 
quantum description of Einstein gravity. It has successfully dealt with tra- 
ditional problems like a macroscopic account of the black hole entropy |56j 
and the construction of non-singular cosmological models [S7], among other 
topics. One of its most notable predictions is the property that area and vol- 
ume operators are quantized in terms of the corresponding powers of £p, thus 
signaling a granular structure of space at short distances j^E] ■ R is precisely 
the possible consequences that such granularity may induce at energy scales 
of the standard model what is to be explored in this work. This question 
is intimately related to the still open problem of the semiclassical limit in 
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LQG: how does one recover the continuous metric description of space-time 
starting from the quantum version of it?. 

2.1 The passage to the new variables 

LQG is formulated in terms of the Ashtekar-Barbero variables that arise as a 
canonical transformation from the usual ADM variables. Here we summarize 
the main steps leading to this choice and closely follow Ref . 55. . 

Let us start from the standard Einstein-Hilbert action for gravity in the 
signature (-, +, +, +) 

S[g^} = ±- K j ( dix ) V C ^%-) (4)i? > K = > ( 8 ) 

where ^'R is the four dimensional Ricci scalar. The ADM variables result 
from a foliation of spacetime in 3D surfaces E of constant parameter t, de- 
scribed by coordinates x a ,a = 1,2,3, together with the parameterization of 
the ten variables in terms of the six components on the induced three- 
metric q a b on E, plus the three components of the shift vector A a and the 
lapse function N. In terms of these variables the four dimensional metric is 

9tt = qcdN c N d - A 2 , g ta = g at = q ac N c , g ab = q ab (9) 

and the action J5J reads 

S[q ab , N a , N] = ±-jdtf (d 3 x) Vdct(q ab ) N ^R + K ab K ab ~(K a a 

(10) 

where all the indices are now lowered or raised by q ab and its inverse q cd . 
Here ^'R is the Ricci scalar of the manifold E. The information about the 
velocities q ab is contained in the extrinsic curvature 

K a b = ^ (qab - C-N qab) = K ba , (11) 

where Cn denotes de Lie derivative along the vector A a . Introducing the 
canonical momenta 

Ijjab = M = I ( R ab _ q ab K c^ ^det(^), (12) 
K Oqab K 

SS 5S 

n = = o, n a = — ^- = o (13) 

SN SN a V ' 

and making a Legendre transformation one finds 

S[q ab , n ab , N a , N] = j-JdtJ (d 3 x) [n ab q ab - N a H a -NH]. 

(14) 
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The above action shows that gravity is a fully constrained theory with zero 
canonical Hamiltonian. The invariance under local diffeomorphisms in S is 
generated by the constraints H a , while the dynamics is generated by the 
Hamiltonian constraint H. These are first class constraints which can be 
written as explicit functions of the canonical variables. The shift and the 
lapse functions turn out to be Lagrange multipliers. 

The symplectic structure corresponding to (|14|l is given by the Poisson 
brackets 

{n ab (x), q cd (y)} = 2^ c S b d / 3 \x,y), (15) 

with the remaining ones being zero. 

The second step towards the construction of the new variables is to intro- 
duce a non-abelian SO(3) w SU(2) formulation of the action (|14|l . To this 
end one rewrites the three-metric in terms of a triad e l a such that 

9a& = e' a 4%, (16) 

where the indices i, j = 1,2, 3, transform under local rotations. 

The canonical variables arising in this step are the densitisized (weight 
+ 1) triad 

E? = l -e abc e l]k e{et (17) 
together with the projected extrinsic curvature 

K a = K ah E)^. (18) 
In terms of them we can rewrite 

det(g ab ) = det^f), q ab II ab = 2K l a Ef , (19) 

and the action is 



S 



Ef, Kl,N a ,N,N l 



= ^J dt J i d3 x) [E1K-N a H a {E,K) 

—NH(E,K) — N i G l (E,K)] , (20) 



The introduction of an extra SO(3) gauge freedom in l|16f) requires the addi- 
tional constraints 

Gi{E,K) = e iik E^K k a: (21) 

which turn out to be just the non-abelian Gauss law. 

The final step to the Ashtekar-Barbero variables starts from the recog- 
nition that there is a natural SO(i) ~ SU(2) connection r* compatible with 
the triad, such that 

+ e VM = °- ( 22 ) 

Moreover, the triad introduced in Ea. 1)16(1 can be expressed in terms of the 
densitized triad as 
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e. = . : £ a bc£ EaEu, e) = ; - K- . (23) 

The final Ashtekar-Barbero connection, which allows to write the constraints 
in a simpler form, is 

4 = r* + 7 ic. (24) 

Here 7 is the Immirzi parameter and the corresponding field strength 
two-form is 

F l = dA* + j Jk A* A A k , F% = e'h F* b . (25) 
The final action results 



Ef, A\, N a , N, N' 1 



«./ / E ? A i~ Na Ha{E,A) 

-NH{E,A)-N i G i {E,A)], (26) 
with the canonical variables satisfying the non-zero Poisson brackets 

{E b j (y),Ai(x)} = K 7 S h a d i j 5^(x,y). (27) 
The explicit form of the constraints are 

d = D a Ef, H b = E? F*, - (1 + j 2 )Kl G u 



H 



(F«-2(l + 7 2 )^ ^) 25? 25). (28) 



The choice 7 = 1 corresponds to the original Ashtekar variables, which de- 
fine a complexified version of Einstein gravity together with a very simple 
structure for the Hamiltonian constraint. The price to be paid is the need 
to incorporate some reality conditions in order to recover the corresponding 
real formulation. 

Following the standard Dirac quantization procedure arising from the 
action together with the Poisson brackets (|2*7|l in the coordinate space 
defined by the connection A\, we promote the canonical variables to operators 
in such a way that 

Ai*[A}= Ai*[A], E?*[A] = -ih K1 -^V[A] (29) 

and impose the operator version of the constraints as null conditions upon 
the wave function. 



2.2 Holonomies and fluxes 

Next we describe the fundamental operator variables employed to formulate 
LQG and which supersede the quantum version of the previous canonical 
variables 1)17(1 and 1)24(1 . The basic support for the operators in LQG are open 
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Fig. 1. Open and closed curves in £ 



or closed curves in S, as depicted in Fig.^ respectively, together with two- 
dimensional surfaces. These allow to define appropriately smeared versions of 
the operators 1|17[) and l|24ll . The SU(2) connection operators A l a (x) are re- 
placed by holonomies (or parallel transport matrices) U along open or closed 
curves in the thee-manifold S, defined as 



The notation is the following: the Pauli matrices a % are the generators of the 
SU (2) Lie algebra. The open curve 7 is parameterized by s' with < s' < s 
and V denotes the standard path-ordered product. The holonomy, which is a 
multiplicative operator, transforms covariantly under SU (2) gauge transfor- 
mations at the end points, but still it is not invariant under diffeomorphisms. 

The canonically conjugated momenta operators Ef are replaced by their 
corresponding fluxes over surfaces S in E 



Here er = (a 1 , a 2 ) are the intrinsic coordinates of the surface and n a is its 
normal vector. The above definitions (|3U[) and i|31fl allow to calculate the 
corresponding commutator, which turns out to be a non-canonical one. 

Quantum states \Prj{A) are represented by functionals of generalized 
connections (SU(2) group elements) defined over graphs in S, which are 
called cylindrical functions and define the kinematical space of the problem. 
As represented in Fig. |3 a graph r = { Vi, . . . , V n ; 71, ... , 7 m } is a set of 
points { Vi, . . . , V n } € S, called the vertices, joined by curves 71, ... , 7™, 
called the edges of the graph. The number of edges attached to a vertex 
is called the valence of the vertex. To each edge of the graph we associate 
a group element U (7, A) labelled by an irreducible representation of SU (2) 
and consider a function / which is a map from the direct product of group 
representations to the complex numbers, so that 

& rj (A) = f(U^n,A), . . . , U(-y m , A)) e Kinematical Space. (32) 

For a given graph r we define the scalar product of two associated functions 
/ and g as 




(30) 




(31) 
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Fig. 2. Piece of a graph F = {Vi, . . . , V„;ji, . . .,7m} 



< &r,f\¥r, B >= / ([[dU^fiUx . . . U m )g{U x . . . U m ), (33) 
J ft 

where we have denoted [/$ :— Ufa, A). Here dUi is the corresponding Haar 
measure for the group SU(2). 

There is still a long way to construct the physical space. To this end one 
has to subsequently impose SU(2) gauge invariance, three dimensional dif- 
feomorphism invariance and finally the Hamiltonian constraint. The first two 
conditions can be explicitly accomplished. The first step is implemented by 
an adequate choice of the functions / in terms of the so-called intertwiners 
at each vertex. These are invariant tensors in SU(2) which map the product 
of representations at each vertex into gauge invariant expressions. In these 
way the so called spin network states are defined and it is possible to prove 
that they form an orthonormal basis for the SU(2) gauge invariant kinemat- 
ical space. The second step is performed by going from loops to knots, thus 
defining equivalent classes of loops under diffeomorphisms. 

The action of the operators (|30|l and l|31|l upon quantum states (|33|l is 
defined by 

(Ui(A, 1 )*)[A] = u£(A, 1 )x*[A], 

Ei(S)U(A,>y) =Y,±MU(A,^)o i U{A, 1 ?»). (34) 
k 

Thus, holonomies act as multiplicative operators, while the action of the 
smeared conjugated momentum depends upon a given edge 7 of the graph 
crossing or not the associated surface S at the points For each intersection 
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Pk the path 7 is separated in two pieces ^ k and 7^* such that 7 = 7^ o 7^* . 
The action at that intersection insert the generator between the holonomies 
corresponding to each of these paths. This is shown schematically in Fig. 
To close this summary we mention the area operator A which is defined in 




Fig. 3. Conditions for the action of the smeared momentum operator 
the following way 

A = lim V JEi(S n )Ei(S n ), (35) 

n— >oo * — » v 

n 

by partitioning the associated surface S. A similar construction is made for 
the volume operator V . These are well defined partial observables (hermi- 
tian operators in the SU(2) gauge invariant kinematical space) with discrete 
eigenvalues and eigenfunctions given by the spin network states. An impor- 
tant property to be used in the following is that only the nodes contribute to 
the action of the volume operator upon states defined in a graph. 

2.3 Elements of Thiemann regularization 

Thiemann has proposed a general regularization scheme that produces a 
consistent mathematical definition for the operators entering in the descrip- 
tion of loop quantum gravity E21 • Such regularized operators act upon 
states which are functions of generalized connections defined over graphs. 
The regularization procedure is based upon a triangulation of space which 
is adapted to each graph. This means that the space surrounding any vertex 
of r is filled with tetrahedra A having only one vertex in common with the 
graph (called the basepoint V(A)) plus segments Ui(A),I = 1,2,3, start- 
ing at V and directed along the edges of the graph. For a vertex of valence 
greater than three one must consider all possible combination of three edges 
to build the tetrahedra. In the regions not including the vertices of r the 
choice of tetrahedra is arbitrary and the results are independent of it. The 
open path along the segment U j is denoted by 7v r ,[/ / • The arc connecting 
the end points of Ui(A) and Uj(A) is denoted by au(A) and the loop 
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A first step in the implementation of the regularization method is to 
express each connection and field strength in terms of an adequate holonomy. 
For example, using the path jv,Ui , of coordinate length e, we can incorporate 
a connection A a (x) into the corresponding holonomy via the expansion 

U(rrv,Ui,A) = l + eU?A a (x) + 0(e 2 ). (36) 

An analogous calculation for the closed triangular path 0^,1/2,1/3 leads to 
the following expression incorporating the curvature 

U(av,u„u a ,A) = l + ~e 2 [/ 2 C/| F ab + 0(e 3 ). (37) 
A second ingredient of the method is the basic identity 61 



{A*,V(R)} PB , V= / (d 3 x)y/det(E ai ), (38) 



y/det{E{x)) 

where R is a region such x € R. In this way the volume V, which will 
be promoted to the operator level once the quantization is performed, is 
introduced in the regularization. 

To fix ideas let us look at the first piece of the smeared Hamiltonian 
constraint 

Hi = I (d 3 x) N(x) / (d 3 x)N(x) e abc Fl b {A a , V}), (39) 



y/det(E) 
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where we have used Ea. (|38fl in the second step. Now, substituting the ex- 
pressions ||5d]) and i|T7|l for A a , F ab € SU(2), respectively, and replacing the 
integral with a sum, we can write 

Hi = lim^e 3 N m e IJK -^Tr( {U{a VmtUlVl ,A)~U- 1 {a v ^ Ul ^ Uj ,A)) 

m 

x U- x { 1Vm , UK ,A) {U( 1Vm ,u K ,A),V(R m )} ) . (40) 

The crucial fact is that the divergent contribution arising from 1/ y/det(E), 
which goes like 1/e 3 , is cancelled by the correspondent factors contained in 
expansion of the holonomies. In this way Eq. 1)40(1 leads to a well defined 
regulated expression. Furthermore, when the correspondent quantities are 
quantized, the volume operator V is naturally introduced. This has the im- 
portant consequence that the action of such operators upon wave functions 
defined on a graph get contributions only from the nodes of the graph. 

Another illustrative example of this regularization is the magnetic sector 
of the electromagnetic Hamiltonian 

l/**8scp*< <41) 

which is translated into the operator expression 

uB _ s J KL MNP v 

~ 2£ 4 O 2 ^ ^ 

p vev(r) v{A)=v{A')=v 

- l) , (42) 

where 

(43) 



Wma = Tr yr k h Ul{A) 
Here we have simplified the notation by introducing 

huj(A) = U(-fv A ,Ur,A), haj K (A) = U (av A ,Uj,U K ,A). (44) 

The underlined quantities (holonomy and connection) refer to the electro- 
magnetic sector, and are to be distinguished from the gravitational ones. 

In the following we will separately consider each matter contribution to 
the Hamiltonian constraint as the corresponding Hamiltonian for the respec- 
tive gravitationally coupled sector. 



3 A kinematical estimation of the semiclassical limit 

The construction of states that live in the physical space of LQG and which 
approximate a given geometry at large distances, while retaining their gran- 
ular structure at Planck length scales is still an open problem. This is usually 
referred to as the semiclassical limit of LQG and it is actively under investi- 
gation 
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3.1 Heuristic characterization of the states 

In order to make some preliminary steps towards the study of the conse- 
quences of space granularity at large scales in the flat space limit, we take 
here an heuristical point of view, starting from the exact and well under- 
stood operator version of LQG. The approximation we introduce consists in 
defining their action upon the semiclassical states through some plausible 
requirements, without having an explicit construction for such states. We 
think of the semiclassical configuration corresponding to a particular mat- 
ter or gauge field, described by operators collectively denoted by F(x), plus 
gravity as given by an ensemble of graphs r, each occurring with probability 
P(r). To each of such graphs we associate a wave function \r, C, F) := \r, S) 
which we assume to be peaked with respect to the classical matter field con- 
figuration F(x), together with a flat gravitational metric and a zero value for 
the gravitational connection at large distances. Not surprisingly, the semi- 
classical approximation gives us simultaneous average information regarding 
a field together with its canonically conjugated momentum. In other words, 
the contribution of the gravitational and matter operators inside the expec- 
tation value is estimated as (21 132 E3| 
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(45) 



and similarly for any product of operators inside the expectation value. As 
an illustration, the gauge field is taken to be the electromagnetic field in the 
following but an analogous estimation will hold for any other gauge or matter 
fields. The parameter T > is a real number and F_ a b denotes the classical 
electromagnetic field strength. We further associate the effective hamiltonian 

H r = (r,S\H r \r,S) (46) 

to each graph. 

The coarse graining scale C » ip of the wave function is such that a 
continuous flat metric approximation is appropriate for distances D much 
larger that C, while the granular structure of space becomes relevant when 
probing distances D smaller that C. That is to say, C is not a scale of quantum 
gravity but rather it is a scale that separates the continuum description of 
space from its discrete characterization that is fully manifest at the Planck 
length £p, which signals the quantum behavior. Summarizing, we expect the 
following behavior 

D >> C >> £p : continuous flat classical geometry, 
l P ~ D « C : manifest discrete structure of space, 
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for the probe scale D. The coarse graining scale C is not provided by our ap- 
proximation and has to be estimated in each particular case. We will explore 
some alternatives in subsection 13.41 



Fig. 5. Different scales for the model of space considered in the text 

As indicated in Fig. [3 we think of space as constructed by adding boxes 
of size C 3 , which center represents a given point x in the continuum limit and 
which contain a large number of vertices of the adapted graphs. The matter 
field F(x), characterized by a de Broglie wave length A, is considered as a 
slowly varying function within each box (A > C) and contributes with its 
classical value at the center of the box, when taking expectation values. The 
requirement A > C guarantees that we can describe the flat space dynamics 
in terms of an effective field theory, using the standard differential calculus of 
the continuum. On the other hand, gravitational variables are rapidly varying 
inside each box. 

The total effective Hamiltonian is defined as an average over the graphs 

r 



3.2 The calculation 

We summarize now the method of calculation 1321 IS'H] . For each graph 
r the effective Hamiltonian is defined by (|46|l . For a given vertex V inside 
the expectation value, one expands each holonomy or flux operator in powers 
of the segments Uj(A)) of the attached tetrahedra A, plus derivatives of 
the gauge or matter fields operators. Schematically, in the case of this 
produces 





(47) 



r 




£ (r,s\(d^d a *...d a *...F piqi (v)) 



. . . x 



ve V(r) v(A)=v 
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x (... d^d^F^V^fvi^iVV^AlA^V) \r,S). 

(48) 

Here T contains gravitational operators (connection and volume operators, 
for example) together with contributions depending on the segments of the 
adapted triangulation in the particular graph. Next, according to Fig.0 space 
is considered to be divided into boxes, each centered at a fixed point x and 
with volume £ 3 ss d 3 x. The choice of boxes is the same for all the graphs 
considered. Each box contains a large number of vertices of the semiclassical 
state (C » £p), but it is considered as infinitesimal in the scale where 
the space can be regarded as continuous. The sum over the vertices in ijl'gl 
is subsequently split as a sum over the vertices in each box, plus a sum 
over the boxes. Also, one assumes that the gauge or matter operators are 
slowly varying within a box (C « A, with A been the correspondent particle 
wavelength), in such a way that for all the vertices inside a given box one 
can write (r, S\ . . . F afc (V) . . . \r, S) — ^F_ ab {x), for example. Here F_ ab is the 
classical electromagnetic field at the center of the box and /i is a dimensionless 
constant which is determined in such a way that the standard classical result 
in the zeroth order approximation is recovered. Applying the procedure just 
described to (|48|l leads to 

#r=E (d ai d a *...d a *...F piqi {x)) ...(... d a ^d a -F pq (x))x 

Box 

J2 ^ E » n+ \r,s\j r f ai ... am ™---^(v,u(A),A ia )\r,s), 

VeBox V(A)=V p 

(49) 

where n + 1 is the total number of factors F rs (x), each of which can contain 
some derivatives. The expectation value of the gravitational contribution is 
supposed to be a rapidly varying function inside each box. Next we consider 
the total effective Hamiltonian (|47|l . which is defined as an average over the 
graphs r, i.e. over the adapted triangulations. This effectively amounts to 
average the remaining expectation values in each box of the sum l|49|l . We 
call this box-average t ai ,,, am pq --- piqi (x) and define it by 

£p(r) E l ? E (r,s\±f ai ... am ^^\r,s) = 

r veBox v(A)=v p 

d 3 xt ai ... am pc >- p ^(x). (50) 

We estimate t ai ... am ^—piQi (gj) by demanding it to be constructed from the 
flat space tensors 5 a b and e a b c , together with the corresponding ones Ti and 
5°, in the case of spinors. In this way we are imposing isotropy and rota- 
tional invariance on our final Hamiltonian in the frame selected by the choice 
of the semiclassical states. This is somewhat analogous to the spontaneously 
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symmetry breaking scheme, where the choice of the vacuum selects a partic- 
ular frame for the symmetry breaking vacuum expectation values. Also, the 
scalings given in (|45|l together with the additional assumptions 

(r,s\...v...\r,s)-^f P , u?-^£p, (5i) 

are used in the above estimation. After replacing the summation over boxes 
by the integral over space, the resulting Hamiltonian has the final form 

H B = J d A x (d ai d a * ...d a K..F pigi (x)] ... ( ...a°»- i a a -'£ M (a:)' 

xt ai ...a m pq - piqi (x)- (52) 

Finally we exhibit some of the required expansions. The basic electromag- 
netic holonomy around the path a JK 

h ajK (A)=e Ja -"< lA) =e I JK >, (53) 

can be written in terms of the magnetic flux $> B through the face bounded 
by such closed path. The expansion of such flux with respect to the magnetic 
field at the vertex is 

$ B (Fu) = U + \{Uj + Wj) d c + ±(U? Uf + Uj Uj + U c j Uj) 

x d c 8 d + ...) ~U?U b je abc B c (V), 

to second order in the derivatives. 

The expansion of the gravitational operator l|43|l is 

Wi la = U a L w la + U a L U b L w lab + UlU b L U c L w labc + 0(s 4 w), (54) 

to third order in the vectors Z7£. The remaining operators in the expansion 
include the connection together with the volume operators, in the form 

w ia = ^[A ia ,W}, Wiab = ~eijk[Aj a , [A kb ,W\], 

Wiabc = -jslAja, [Ajb, [Ac, W}]]. (55) 



3.3 The results 



The application of the method described in the previous section leads to the 
following effective Hamiltonians. In the electromagnetic case we obtain [32] 
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(56) 



The above result has been extended to the Yang-Mills case in Ref.|S3]- The 
case where only 9s in non zero corresponds to the Gambini-Pullin effective 
Hamiltonian [M)\ . 

In the case of Majorana fermions with mass m we have OS] 
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r d d d i{x) 
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r a a a e(x) + c.c. 



(57) 



where ?r(x) = i£*(x) is the canonically conjugated momentum of the two- 
component spin 1/2 field £(x). The notation is r a = —(i/2)a a , where cr° are 
the standard Pauli matrices. 

Some comments regarding the procedure are now in order, (i) The di- 
mensionless numerical coefficients 9a, ka appearing in the above effective 
Hamiltonians remain arbitrary in the procedure, but are independent of any 
parameter of the model. To predict them will require an exact knowledge of 
the states used to calculate the corresponding average values, together with 
the action of the basic operators upon them. In this respect we have obtained 
a systematic parameterization of the possible modified Hamiltonians in terms 
of higher derivative operators, where the dependence upon the scales of the 
problem (£p,£,T) has been explicitly determined. In particular, these coef- 
ficients could turn out to be zero, leading to no dynamical corrections, (ii) 
A main drawback of the method is that it does not incorporate properly the 
dynamics via the Hamiltonian constraint of quantum gravity. An improved 
semiclassical approximation taking care of the many issues assumed or left 
over in our approach is certainly needed, (iii) There is also the question of 
interpreting the results H56JI and Ij57(l in relation to their transformation prop- 
erties under standard active (particle)global Lorentz transformations. There 
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are various possibilities. These Hamiltonians have been calculated in a spe- 
cific reference frame where isotropy under rotations is maintained. This could 
signal the presence of an absolute reference frame, a reborn version of the 
ether. This hypothesis will be further explored in the following and has re- 
ceived a lot of attention in the framework of different models in the literature. 
On the other hand one could argue that this specific frame has been selected 
by the choice of the semiclassical states that describe the fiat space limit, 
in analogous way to the spontaneously symmetry breaking procedure. The 
transformation properties under this assumption have not been studied. The 
third alternative is that only a complete calculation of the corrections to flat 
space dynamics will tell us which is the proper generalization, if any, of the 
global Lorentz transformations, thus yielding a deformed or extended version 
for preserving the equivalence of inertial frames and leading to a modified 
relativity principle, (iv) Finally we must say that the above results have pro- 
vided some motivation for phenomenological theories which have explored the 
consequences of such modifications and have shown that either experimen- 
tal or astrophysical observations, even with the actual level of sensitivity 
can set rather stringent bounds upon the correction parameters. In other 
words, Planck scale sensitivity has been already attained in a rather broad 
set of observations. Conversely, this would imply severe constraints upon the 
fundamental theory once the semiclassical limit is correctly performed. 

3.4 The parameters C and T 

In order to produce numerical estimations for some of the effects arising from 
the previously obtained modifications to flat space dynamics, we must further 
fix the value of the scales C and T. 

Let us recall that £ is a scale indicating the onset of the distance from 
where the non perturbative states in the loop representation can be approx- 
imated by the classical flat metric. The propagating particle is character- 
ized by energies which probe distances of the order of the De Broglie wave 
length A. As mentioned previously, just to be consistent with a description in 
terms of classical continuous equations it is necessary to require that C < A. 
Two distinguished cases for choosing C arise: (i) the mobile scale, where 
we take the marginal choice C — A in each situation and (ii) the universal 
scale, which has been introduced in the discussion of the GZK anomaly |5U). 
The consideration of the different reactions involved produces the preferred 
range 4.6 x 1(T 8 Gey- 1 > C > 8.3 x lO^GeV" 1 . A recent study of the 
gravitational Cerenkov effect together with neutrino oscillations [HE] yields a 
universal scale evaluation which is consistent with the former. 

Ranges for T have been estimated considering the observation that atmo- 
spheric neutrino oscillations at average energies of the order 10 -2 — 10 2 GeV 
are dominated by the corresponding mass differences via the oscillation length 
L m |fifi| . This means that additional contributions to the oscillation length, 
in particular the quantum gravity correction Lqg, should satisfy Lqq > L m . 
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This has been used to set lower bounds upon T. Within the proposed two 
methods for estimating the scale C of the process we obtain: (i) T > 0.15 
when C is considered as a mobile scale and (ii) T > 1.2 when the scale C 
takes the universal value C « 10 -8 GeV~ x [33]. These results are summarized 
in Table E 



Table 2. Ranges for the parameters £ and T 
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4 Phenomenological aspects 

In this section we discuss four phenomenological applications of models pre- 
senting Planck scale dynamical corrections. In subsection 14. II we summarize 
the derivation of stringent bounds upon combinations of some parameters ap- 
pearing in the effective fermionic Hamiltonian (|57fl , using existing data from 
clock-comparison experiments. Subsection l4.2l is devoted to the calculation of 
radiative corrections arising in LIV effects with preferred frames associated 
to space granularity. In order to make the calculation simple we consider a 
Yukawa model for the interactions. Subsection 14. 31 contains the discussion of 
synchrotron radiation in Planck scale modified electrodynamics, emphasizing 
the model of Ref. which is phenomenologically formulated as a theory pa- 
rameterizing LIV with dimension five operators. Finally, a unified description 
of radiation in Planck scale modified electrodynamics incorporating different 
models is presented in subsection 14.41 

4.1 Bounds on the fermionic sector parameters from 
clock-comparison experiments 

The Hamiltonian (|57|l for two components fermions was obtained under the 
assumption of flat space isotropy and it was assumed to account for the 
dynamics in a preferred reference frame, identified as the one in which the 
Cosmic Microwave Background looks isotropic. The earth velocity w with 
respect to that frame has already been determined to be w/c~ 1.23 x 10~ 3 
by COBE [S7]. Thus, in the earth reference frame one expects the appear- 
ance of signals indicating minute violations of space isotropy encoded in the 
diurnal variation of the w-dependent terms appearing in the transformed 
Hamiltonian or Lagrangian |fi8| . On the other hand, many high precision 
experimental test of such variations, using atomic and nuclear systems for 
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example, have been already reported in the literature ^3 E3 d] an d 
already analyzed in terms of the SME. Amazingly enough such precision is 
already adequate to set very stringent bounds on some of the parameters 
arising from the quantum gravity corrections. 

We have considered the case of non-relativistic Dirac particles and ob- 
tained corrections which involve the coupling of the spin to the CMB ve- 
locity, together with a quadrupolar anisotropy of the inertial mass |68| . The 
calculation was made with the choices T = and C = 1/M, where M is 
the rest mass of the fermion. Here it is important to emphasize again that 
C >> lp is not the scale of quantum gravity but rather the lower limit of 
distances from which a description of space as a continuum is already valid. 
In this sense we are taking the upper limit settled by the de Broglie wave 
length Xd — 1 /M of the nucleon, which allows us to make sense of the corre- 
sponding Dirac equation. As we will show later, a lower limit for C will only 
make the bounds upon the LIV parameters more stringent. 

Keeping only terms linear in tp, the equations of motion for a Majorana 
fermion of mass m described by the two-component spinor £, arising from 
the Hamiltonian 157J1 are given by 

[id/dt - iAa • V + K/2] £ - m{a - i(3cr ■ V)x = 0, 

[id/dt + iAer-V- K/2] x - m(a - i/3cr ■ V)£ = 0, (58) 

where 

A = (l+G 1 me p ), a= {l + 3 m£ P ), 

K = mOimip, f3 = 2 £ P . (59) 

The notation in Eqs. (|58(l is x — ~*02£* ana * ^ ^ s a direct matter to verify 
the consistency between them. These equations can be readily extended to 
the Dirac case by considering \ an d £ as independent spinors unified in 
\p T = (£ T ; x T ), with the result 

+ e im £ P i-y-V- y757° - m (a - i0 2 £p £-V)J& = 0. (60) 

Here we have used the representation in which 75 is diagonal and the spin op- 
erator is S k = (i/2)ekiml l l m , with standard particle physicist conventions. 
The normalization has been chosen so that in the limit (m£p) — > we re- 
cover the standard massive Dirac equation. The term m (1 + ©3 m£p) can be 
interpreted as a renormalization of the mass whose physical value is taken to 
be M = m(l + ©3 mlp). After these modifications we can write an effective 
Lagrangian describing the time evolution as seen in the CMB frame. In order 
to obtain the dynamics in the laboratory frame we implement an observer 
(passive) Lorentz transformation in the former Lagrangian and rewrite it in 
a covariant looking form by introducing explicitly the CMB frame four ve- 
locity =7(1, w/c) arising from the boosted rest frame velocity (1,0). 
The result is 



Dynamical corrections from space granularity 23 



+ j(0 2 M£ P )'Pe fwa pW^ 1 a d^ - j(9 4 Ml P )MW^ 5 ^& + h.c, 

(61) 

where is an external non-dynamical field which is not to be varied in the 
corresponding action. It is interesting to remark that the above Lagrangian 
provides a specific realization of the general form considered in the SME. 
According to Ref. the identifications are 
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(62) 



The above reference provides also the non-relativistic limit of the Hamiltonian 
corresponding to the general SME modified fermion Lagrangian. In our case, 
up to first order in lp and up to order (w)/c 2 , the identifications (|62[1 yield 



H = Mc 2 (l + 6>i M£ P (w/c) 2 ) + OxMli 
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(63) 



where s l = a % /2. 

The above effective Hamiltonian has been employed in the description 
of the valence nucleon responsible for the transitions measured in clock- 
comparison experiments using pairs of nuclei like ( 21 Ne, 3 He) 70, and 
( 129 Xe, 3 He) [7T], for example. In (|63|l we have not written the terms lin- 
ear in the momentum since they average to zero in the nuclear bound state 
situation. Here g is the standard gyromagnetic factor, and Qp is the momen- 
tum quadrupole tensor with components Qpij = PiPj — l/3p 2 Sij. The terms 
in the second square bracket of the LHS of iffl^ represent a coupling of the 
spin to the velocity of the privileged reference frame. The first term inside the 
bracket has been measured with high accuracy and an upper bound for the 
coefficient has been found. The second term in the same bracket is a small 
anisotropy contribution and can be neglected. Thus we find the correction 
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(64) 



The first square bracket in the LHS of l|63() represents an anisotropy of the 
inertial mass and has been bounded in Hughes-Drever like experiments. With 
the approximation Q p = —5 / '3 < p 2 > Q / R 2 for the momentum quadrupole 
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moment, with Q being the electric quadrupole moment and R the nuclear 
radius, we obtain 

W„ («)(=)•«(»»), (65) 

for the quadrupole mass perturbation, where 9 is the angle between the quan- 
tization axis and w. Using < p 2 /2M >~ 40 MeV for the energy of a nucleon 
in the outer shell of a typical heavy nucleus, together with the experimental 
bounds of references |7()l ITT] we find [HE] 

I ©2 + ^G>4 |< 2 x 1CT 9 , | 6>i |< 3 x 1CT 5 . (66) 

From Eqs. and ljfi5]) we realize that only products of the type (LIV pa- 
rameters) x l/£ k (LIVparameters) x M are bounded by the experiments. 
In this way choosing a smaller C, i.e. a larger M, will just decrease the upper 
bound for the LIV parameters. Thus the choice C = Xd = 1/M produces the 
weakest bound in this analysis. 

4.2 Radiative corrections in preferred frames modelling space 
granularity 

So far corrections to the dynamics have arisen only in the non-interacting 
theory, through factors of the type (E/Mp) , which are directly relevant at 
unaccessible energies E ss Mp. A possibility of probing such high energies 
is through the inclusion of radiative corrections (particle's self energies, for 
example ) because the internal momenta are integrated up to the maximum 
allowed in a given reference frame. The standard folklore with respect to 
any new physics entering at high scales (here Planck scale) is that it has 
negligible effects on the leading-order low-energy physics (here free particle 
corrections). Contrary to this belief, we show in this subsection that mod- 
elling space granularity via the introduction of a physical cutoff, which defines 
a preferred reference frame, leads to unsuppressed dimension four LIV con- 
tributions |74l ITB*] . Results consistent with these have been obtained in Ref. 
d 

To this end we consider the calculation of one-loop radiative corrections 
in the Yukawa model 

L = 1/2 (c^) 2 - 1/2 /i 2 </> 2 + ^(ij^dfj, — m)ip + gcfnptjj + (LVT). (67) 

Here LVT refers to the highly suppressed zeroth order Lorentz violating terms 
that take into account the previously discussed free particle dynamical mod- 
ifications. We model the space granularity by introducing a physical cut-off 
A in such a way that the magnitude of the three-momentum in any loop is 
bounded by this quantity. The parameter A defines the onset of the scale 
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at which the granularity of space becomes manifest. A convenient way of 
incorporating this requirement is to introduce the physical cutoff function 

f(\k\/A), /(0) = 1, /(oo) = 0, Anl/C, (68) 

which suppresses the internal momenta having \k\ > A. One can visualize 
this function as a smoothed theta function with the correspondent jump at 
|fe| w A. Our choice is that the cutoff function depends on spatial momentum, 
but not of energy, and that it is rotationally invariant. 

The fermion self energy \75\ 

Tests of Lorentz invariance typically deal with the relation between energy 
and 3-momentum of an isolated particle, which is given by the position of 
the pole of the particle's full propagator. We will calculate the effect of loop 
corrections on this relation. It is convenient, as usual, to write the full fermion 
propagator as 

Sip) = -z v( y (69) 

where S% (p) is the standard fermion self energy. Here we are neglecting the 
zeroth order corrections to the dispersion relations for the free particles. 

As a first manner to model the validity of a continuous description of 
space only up to some short distance 1/ A » £, we choose to cutoff only the 
free scalar propagator, but not the free fermion propagator, which is sufficient 
to cutoff the UV divergence of the one-loop fermion self-energy while giving 
a maximally simple calculation. In this way, the one-loop approximation to 
the fermion self-energy is 

. 2 f d 4 fc f(\k\/A) lp (p + ky + m 

where we have set fi = for simplicity. We wish to investigate the self-energy 
when the momentum p^ and the mass m are much less than the parameter 
A. If the integral were convergent as the parameter is removed (A — * oo), 
then the self-energy would equal its value at A = oo plus corrections of order 
m/A and p/A. However, it is easy to see that the integral is divergent in this 
limit. To analyze its behavior, we observe that differentiating it twice with 
respect to p gives a finite integral (when A — > oo). So we extract the zeroth 
order and the linear terms in the momentum expansion of S 2 about p = 

27 2 (p) = Z 2 (p = Q)+p d p oS 2 (p)\ p=0 +p l d p >Z 2 (p)\ p = + 2{p). (71) 

Now S(j)) is 0(p 2 ) as p — > and it is finite when the cutoff is removed. There- 
fore it only gives power-suppressed Lorentz violations. The unsuppressed 
Lorentz violations, if any, can only arise from the first three terms. We find 
that the terms S%ip = 0), d p aS<zip)\ p =a and d p iS2(p)\p=o are proportional 
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to the Dirac matrices 1, 7,, 7 , exactly as in the Lorentz- invariant case; this 
follows from the discrete symmetries of our model, rather than full Lorentz 
invariance. Hence we can write 

S 2 {p) = Am+p Q lQ B + P l ^C + E(p, A = 00) + 0(p/A,m/A) 

= Am+ p^C + p°7o(B - C) + S{p, A = 00) + 0{p/A, m/A), 

(72) 

where A, B and C are numerical-valued functions of the parameters of the 
problem. The only leading-power Lorentz violation is caused by the difference 
between the B and C, which we calculate below and find it to be nonzero 
and unsupprcssed. 

In fact, the constant term is 



7 P fc p + m 



(k 2 — m 2 + ie) (k 2 + ie) ' 



(73) 



Since the physical cutoff depends only on |fc|, the modified theory is invari- 
ant under reversal of any component of k^ and the term proportional to 7 P 
vanishes after integration. Hence 
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where E, 



VP 



In the second line, we have performed the integrals 



over k° and the angle of k. The corresponding term in S 2 is proportional to 
the unit Dirac matrix, so that it is equivalent to a fermion mass term. Note 
that when A — > 00 there is a logarithmic divergence, which can be removed 
by the usual mass renormalization. Since the term is proportional to the unit 
Dirac matrix, it does not violate Lorentz invariance. 

The situation for the derivative terms is different. For the time-like deriva- 
tive we have 



d p o£ 2 (p)\ P =o = 
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(75) 



Invariance under reversal of components removes all but the term propor- 
tional to 70, so that 
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Similarly, from the space-like derivative we get 



C 
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(77) 



As expected, both B and C have logarithmic divergences as A — > oo. 

We are interested in the difference, which gives the low-energy violation 
of Lorentz invariance 
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(78) 
(79) 



Suppose the integral (|78|l were convergent in the continuum limit. Then it 
would be zero in this limit, as can be seen by Wick rotation of an integral of 
^k 2 +k° times an arbitrary Lorentz-invariant function F(k 2 ). This argument 
depends on the integral being absolutely convergent, i.e., on the integral of 
the absolute value of the integrand being finite. In contrast, the corresponding 
integral for B — C is logarithmically divergent. So the leading term for the 
integral, obtained when A — > oo, depends on how this cutoff is removed. In 
our case this is dictated by the physical modelling of the space granularity via 
the function l|68fl , as expressed in Eq. I|79|l . This choice implies that the factor 
in square brackets behaves like m 2 /l 3 when l/m — > oo, so that the integral is 
finite when the cutoff is removed. Besides, we can write this square bracket 
as indicated in the second line of the following equation, which leads to a 
direct calculation of the leading term 
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(80) 



The nonzero value of B — C shows that the one-loop fermion self-energy 
introduces a Lorentz- violating term at small momenta even when the cut-off 
is made very large: Lorentz violation is suppressed not by a power of m/A, 
but only by a factor of the coupling. A convenient interpretation of this result 
uses the language of effective field theories. Observe that the Lorentz violation 
caused by the B — C term is equivalent at the one- loop level to adding a term 



(81) 
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to the Lagrangian. That is, we could obtain the same results in the A S> m 
region by dropping the B — C term in the self-energy (|72|) and then adding 
this extra term to the Lagrangian. In this extra term £ = g 2 /487r 2 , is the 
coefficient calculated above, and = (1,0) covariantly specifies the rest 
frame of the cutoff. Summarizing, the fermion self-energy contribution has 
the effect of feeding back terms of large unnatural size into the Lagrangian, 
thus producing a strong fine tuning problem when higher corrections are 
taken into account. 

The boson self energy 

Similar results can be obtained from the boson self energy. Since the internal 
lines are both fermionic in this case, a more symmetrical way of incorpo- 
rating the physical cutoff is required. This is accomplished by extending the 
fermionic propagators to 

. ami a) (82) 



jk — m + ie jk — m + B. (k) + ie 

The function S'(fe) takes into account the modified dispersion relations of the 
theory to zeroth order and its contribution will not be taken into account in 
the following. 

Consider only the physical momentum cutoff function and calculate the 
boson self-energy II(E,p, A), defined in terms of the full propagator A(E,p) 

A(E > P) p 2 - J? -n(E,p, A)-ie (83) 
The one loop calculation is 

,2 



n{E,p,A) = -^- i J d 4 kf(\k\/A)f(\k + p\/A) 
Tr[(7 ■ A; + m) (7 • (k + p) + m)] 



(84) 



(fc 2 — m 2 + ie) [(k + p) 2 — m 2 + ie] 
The result can be presented as an expansion in even powers of the momenta 
n{E,p,A) = A + p 2 B + mc p»p v W ll W v + n{E,p,A), (85) 

where TI is convergent when the regulator is removed so that it contributes 
only with Lorentz violating terms which are suppressed by powers of p 2 / A 2 . 
In the frame where the cutoff A is defined, W M = (1, 0), we obtain 

1 / ( d 2 n{ P ) \ ( d 2 n{ P ) \ \ 

%c = HlwJ F0 + lwjJ' (86) 

from 1)85(1 . The calculation of the required pieces starts from l(84|l and after 
some algebra produces 
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(87) 
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to the argument of /. As usual we have Ef. 



f(k/A) and the derivatives, denoted by primes, are with respect 

2 Each expression 



\/k 2 + m 2 



H87|) and (|88|l is logarithmically divergent, but their combination (|86|l is fi- 
nite. Thus we can estimate the leading contribution to tjrc by setting / = 1 
inside the integrals with no derivatives with respect to /. The terms contain- 
ing derivatives are handled by assuming that the contributions of /', /" are 
sharply localized in the region x = k/A m 1. This implies that we neglect 
the mass contributions m 2 << k 2 m A 2 in the corresponding integrals. Also, 
we set equal to zero the resulting boundary terms which include /' in the 
required integrations by parts. The final result is 
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(89) 



which is estimated in the range t]rc > 10 -3 using standard model couplings. 
On the other hand, rjnc can be interpreted as a correction 6c to the boson 
(photon) speed, for which extremely tight bounds exists 
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The two examples presented here show that radiative corrections in pre- 
ferred frames associated to space granularity, modelled by the physical cutoff 
function 1)68(1 ■ induce LIV contributions which are enormously increased with 
respect to the already established bounds. That is to say, they are not sup- 
pressed by factors 1 /M as expected initially. In this way, a naturalness prob- 
lem arises and a fine tuning problem appears when considering higher order 
radiative corrections. One possibility to guarantee the stability of dimension 
three and four LIV contributions is via the introduction of a custodial sym- 
metry. This option has been explored using supersymmetry in Ref . |7fi| . 

An alternative way of dealing with radiative corrections in LIV processes 
can be found in Ref. [77] • 



4.3 Radiation in Lorentz violating electrodynamics 

The observation of 100 MeV synchrotron radiation from the Crab Nebula 
has recently been used to impose stringent limits upon the parameters de- 
scribing a modified electrodynamics embodied in Maxwell equations, which 
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together with the corresponding coupled equations for the charges, get cor- 
rection terms which are linear in the Planck length ^H]. Such bounds were 
based on a set of very reasonable assumptions on how some of the stan- 
dard results of synchrotron radiation extend to the Lorentz non-invariant 
situation. This caused some controversy in the literature |78|. Moreover, an 
assessment of such assumptions requires the introduction of specific dynam- 
ical models. One of them is the Myers-Pospelov (MP) effective theory at the 
classical level, which parameterizes LIV using dimension five operators Ref. 
|14|. A detailed description of synchrotron radiation in this model has been 
presented in Refs.|S2], which we review in this subsection and that has moti- 
vated the general point of view summarized in subsection 14.41 given in more 
detail in Ref.jHH]. 

The study of radiation in LIV electrodynamics constitutes an interest- 
ing problem on its own whose resolution will subsequently allow the use 
of additional observational information to put bounds upon the correction 
parameters. For example we have in mind polarization measurements from 
cosmological sources. The case of gamma ray bursts has recently become 
increasingly relevant [701, although it is still at a controversial stage |80) . 

Our calculation of synchrotron radiation in modified electrodynamics rest 
heavily on previous work reported in Refs. [5Tl 1521 15511541 I55| . 

A partial list of previous studies in electrodynamics incorporating LIV 
via dimension three and four operators is given in Ref. 86 

Synchrotron radiation in the Myers-Pospelov model 

This model parameterizes LIV using dimension five operators both in the 
matter and electromagnetic sectors. There is also a preferred frame four 
velocity V^, which is not a dynamical field. As usual the model exhibits 
passive (observer) Lorentz covariance, which means that the fields and the 
four-velocity transform as tensors when going from one observer frame 
to another. On the other hand, in each frame we violate active (particle) 
Lorentz transformations; that is to say we have the non-dynamical physical 
field in the action, in analogy to the physics going on in the presence 
of an external magnetic field which violates active rotational invariance, for 
example. 

The charge sector 

The dynamics of a classical charged particle of mass /i can be obtained from 
the action for a scalar charged field. In this case the Myers-Pospelov action 
is 



with the notation V ■ d — d^. In momentum space, where we write 
(f(x) = tpo expi(p t — p • x) and in the reference frame where V a — (1, 0), 
the modified dispersion relation becomes 




(90) 
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r/lp 



fj, 



(91) 



To make contact with the results in Ref. 15 it is necessary to rewrite 
fj = —r]/M, rj < 0, where rj is a dimensionless constant and M is a mass 
scale characterizing the Lorentz symmetry breaking, which is usually, but 
not necessarily, identified with the Plank mass. The equation (|91|l is an exact 
relation in fj. From here we obtain the Hamiltonian for a massive particle to 
second order in fj 
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H=(p 2 



2\l/2 



+ P 2 )+^ 3 (p 2 + f 2 ) S/2 + Ofi 3 ). (92) 



In the following we consider the interaction of a particle having charge q with 
a static constant magnetic field B ~ V x A. The standard minimal coupling 
produces the Hamiltonian 
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(93) 



to order 0(fj 3 ). Here c = 3 x 10 10 cm s _1 denotes the uncorrected velocity of 
light in vacuum. In the following we set c = 1. Observe that the dispersion 
relation (|91fl provides the exact inversion 



(p - qA) 2 = (1 + f,E) E 2 



(i 



(94) 



with E being the energy of the particle. The Hamilton equations arising from 
l|9"3*|) yield the acceleration 



^=^[ l -lflE+- A fE z \ (vxB 



(95) 



As in the usual case, this means that the magnitude |v| of the particle velocity 
is constant and that the projection of the trajectory in a plane perpendicular 
to B is a circular orbit with a Larmor frequency 
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(96) 



In general the motion is an helix with pitch angle (the angle between the 
velocity and the magnetic field) a. We restrict ourselves to the motion in the 
plane perpendicular to the magnetic field, i.e. a = tt/2. The solution to the 
equations of motion can be written as 



v(t)-- 



— cos woe, — sin wot, 
Uq uj 



(97) 
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where we emphasize that we are using the standard definition (3 = |v|/c. The 
velocity is 

v(t)=(-/3smu t, (3cosuj t, 0). (98) 

The equation (|97|l identifies R — P/ojq as the Larmor radius of the orbit. The 
Lorentz factor 7 is given by 



7 



2 ' E 2 



pi 1 4 ^ 



(99) 



where the range of energies to be considered is such that << 1, fjE << 1. 

According to the preceding analysis, the current for a charged particle 
moving in the circular motion is 

j(i,r)=^ 3 (r-r(t))v(i), (100) 

where r(t) and v(t) are given in Eqs. (|97|l and l|98() respectively. 

The electromagnetic sector 

The corresponding action is ^2] 



^photon I d X 



(101) 

We choose to work in the rest frame V 11 — (1, 0) where the modified Maxwell 
equations are (c = 1) 

V-B = 0, VxE+— = 0, 

at 

V • E = 47rp, -— + V x B + £<9 (-V x E + d B)) = iirj. 

(102) 

For future purposes it is convenient to define £ = £/M, where £ is a dimen- 
sionless parameter. Introducing the standard potential fields = (cj), A) in 
the Coulomb gauge we have 



= - 47r ^ ( 103 ) 



d 2 A , ~ d 2 A ( 1 



V^A + 2(Vx — = 4d J-V^V • j = 4^ j T , (104) 



<9i 2 s at 2 

where the electric and magnetic fields reduce to 



BA 

E = — — , B = VxA, (105) 



in the radiation approximation. 



Dynamical corrections from space granularity 33 

The energy momentum tensor for this modified electrodynamics is 
given by 

T° ^ 1 L(E 2 + B 2 )-i F E.M, (106) 
S =3 LExB-iExf, (107) 

which are exact expressions in £ and satisfy the usual conservation equation 
outside the sources. 

To solve the equation of motion for A, Eq. (|104fl . we can go to the mo- 
mentum space with the convention 

F(t,r)= f -^ e — t+ik r F( W ,k). (108) 
J (27r) 

The different types of Fourier transforms are specified by the corresponding 
arguments. For example, if F(t, r) denotes the function in space-time, F(u>, r) 
denotes the Fourier transformed function to frequency space, while F(u>,\l) 
denotes the Fourier transformed function to frequency and momentum spaces. 
In this way Eq. (|104Jl reduces to 

-cj 2 + fc 2 -2^w 2 kx) A(w,k) =47rj T (w,k). (109) 

This equation can be diagonalized using the circular polarization basis, with 
A = ±1, giving 



-U) 



2 



k 2 =F2fw 2 fc) A ± (w,k) = 4vrj^(w,k). (110) 



The components Cf of any vector C in the polarization basis associated to 
the direction k are 



fe = g ~ kikk + Xieijkkj^J . (Ill) 



(112) 



The simplest way to proceed is by introducing the total retarded Green func 
tion 

and to calculate 

/r/ 3 k / 
— 3e lk -( M )[G re *(^,k)] lfc . (113) 
(27r) 

After some rearrangements the final integration over dk can be performed by 
the method of residues in the complex plane. Finally one can identify back 
the polarization components of the total Green function as 
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^ U ' r ^ = SR^ eta>{,M ' (U4) 

where R = |r — r'|. Here we have introduced the polarization-dependent 
refraction index n\{z) 



(z) = \/l + z 2 + Xz, z = £w. (115) 



In this way, the fields A A in Eq. (|110|l have well defined phase velocities 
v\ = l/n\(z) and this situation can be described as the propagation of 
photons in a dispersive birefringent medium. 

The Green functions (|114fl determine the corresponding potentials with 
the standard replacements 1/R ~ 1 / 1 r | = 1/r in the denominator, together 
with the following expansion of the phase n(Xz)ujR 



nx(z)u |r — r'| ~ lot 



- n ■ r 1/r 



A£w - A£w- 

r 2 



(116) 



where n = r/r is the direction of observation. Notice that we are interested 
in the radiation approximation of the phase l|116fl . which means that the 
subdominant terms of order (r'/r) 2 or higher are neglected. Consistency de- 
mands that the terms proportional to the LIV parameter £ are larger than 
the neglected one in order to properly include them in this phase. Our general 
results are presented in this full far-field approximation. 
Using H114(l we finally get 

A>,r) = - " A(Z) e"^>""j>,k A ) (117) 
r VI + z 2 

in the radiation approximation. The fields A + (w, r) and A~(u>, r) correspond 
to right and left circular polarization respectively. Let us emphasize that the 
momenta 

k x =n x {z)u;h (118) 

in Eq. I|117|) are fixed in terms of the frequency and the direction of ob- 
servation. The full vector potential is given by the superposition A(u>, r) = 
A + (iv 7 r) + A~(u>, r). Hence the electric and magnetic fields are 

B( Wj r) = \^=^ £ An 2 (z)e m ^j>,k A ), 

E(w,r) = --7=^=5 £ n x (z)e m ^ r i x (Lu,k x )- (119) 

Note that, contrary to the standard case where n x E cx B, here we have 

n x E(w, r) = J_ [B(w, r) + izE(lu, r)] . (120) 
VI + 
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The angular distribution of the power spectrum is defined as 

d 2 P{T) 



dtudft 



(121) 



where P(T) is the radiated power at time T into the solid angle dfl. We can 
compute the total energy emitted in terms of the Poynting vector (|107|) 

/oc 
dt n- S(t,r) r 2 dQ. (122) 
-oo 

This last expression can be rewritten introducing the Fourier transform of 
the Poynting vector, 

POO P J2 771 POO 1 

E= du df2 =/ — [n-S(w,r) + n-S(-w,r)]r 2 dJ?, 

Jo J dilduj J 2n 

(123) 

and allows us to obtain the angular distribution of the energy spectrum 

''' E ' [n-S(o;,r) + n-S(-a;,r)], (124) 



dflduo 2n 

from where the angular distribution of the power spectrum can be identified 

as 

d 2 E f +OD d 2 P{T) 

dflduj J-oo dwdQ ^ 

Writing the Poynting vector I|1U7|) in terms of the polarized potentials we 
obtain the intermediate result 

d 2 E : 



Vl + z 2 [A_(-w,r) • A+(w,r) + A_(w,r) • A+(-w,r) 



dJ2du; 47T 2 

(126) 

Next we express the products A T (— w,r) • A±(w,r) in terms of the current 
j(u>, k) via the relation (|117|) . Using the properties of the fields in the polar- 
ization basis, together with the general relation j k {— ui,— k) = j^(cj,k), we 
obtain 

d 2 E 1 u> 2 r o / s .* / , \ r,+ ■ i i \ 

+n 2 _(z) #(w,k_) P fe - > (w,k_) ] . (127) 

In order to identify the angular distribution of the power spectrum we need to 
introduce the time dependence via the corresponding inverse Fourier trans- 
form. Each contribution in Eq. (|127|) is of the type 

/+oo 
dtdt'e-^-^jl^k) X kr j r (t',k). 
-oc 

(128) 
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Changing to new time variables r = t — t' and T = (t + 1') /2 we get 



/-t-oo poo 
dT / dre- iuT j* k (T + r/2, k) X kr j r (T - r/2, k) . (129) 
-oo J — oo 

Inserting this last relation in Eq. (|12T|> and comparing with Eq. Ijl25|) we 
obtain the final expression for the angular distribution of the radiated power 
spectrum 

d 2 P(T) 1 uo 2 



dwdfl Air 2 ^1 + z 

]T nl(z)f k (T + r/2, k A ) j r (T - r/2, k A ) , (130) 

A 

as the sum of the contributions of both circular polarizations. 
Synchrotron radiation 

This corresponds to the choice l|98f) for the velocity of the current 11 UUP in the 
general expression (|130|) . Following the method of Ref . the time-averaged 
angular distribution of the radiated power spectrum is 

' A=+,- m=0 



dP m , x uf n q 2 l 



[\pn x (z m )J^(W Xm ) + cote J m {W Xm )f 



da 4tt yrr^ 

(131) 

written as a sum over the contribution of the harmonics u> m . Here, the average 
() is taken with respect to the macroscopic time T and J m , denote the 
Bessel functions and their derivatives respectively. The argument of the Bessel 
functions is W\ m = mn\{z m )f5 sin 9. 

We also have calculated the total averaged and integrated power radiated 
into the m-th harmonic 

Pm = u h r=! ^ [ J 2m( 2m n x (z m )(3) 



R yl 

1-I3 2 n 2 (z m ) r 2mn ^)P 



2p 2 n 2 (z m ) J 



dx J2m(x) 



(132) 



which clearly indicates the contribution of each polarization P\ m . The above 
result is exact in z m and the parity-violating contribution has vanished after 
the angular integration. 

In Table |3 we present a rough estimation of the relevant parameters asso- 
ciated with some observed cosmological objects. There r [l.y.] is the distance 
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Table 3. Data of some relevant astrophysical objects 



Object 


r(l.y.) 


7 


B (Gauss) 


u c (GeV) 


uo{GeV) 


m 


m/7 


CRAB 


10 4 


10 9 


10~ 3 


w- 1 


10" 


-30 


10 29 


10 20 


(Mfcn 501) p 


10 s 


10 11 


10 2 


10 4 


10" 


■29 


10 33 


10 22 


(Mkn 501) e 


10 8 


10 11 


io- 1 


10 4 


10" 


■29 


10 33 


10 22 


GRB 021006 


10 10 


10 6 


10 4 


10 -3 


10" 


-18 


10 15 


10 10 



of the object to the earth, 7 is the Lorentz factor of the charged particles at 
the cutoff frequency, B [Gauss] is the average magnetic field producing the 
synchrotron radiation, u c [GeV] is the cut-off frequency and u>o [GeV] is the 
Larmor frequency. In all cases the cut-off frequency lj c has been estimated 
from the radiation spectrum fitted by a synchrotron model in the correspond- 
ing reference. This, together with the magnetic field B allows us to estimate 
the Lorentz factor 



o / ujJGeV] M , . 

7 = 2.36xl0 8 W-^ L 133 

W B[Gauss\ m e 

where M is the mass of the charged particle. From the above we further 
obtain the Larmor frequency 

oj Q [GeV] = 0.6 x 1CT 17 (^f\ B[Gauss\. (134) 

In the case of CRAB Nebula we adopt the typical values given in Ref. |87]. For 
Mkn 501 we consider two possible models for synchrotron radiation where the 
emitter particles are either protons [EE] or electrons [EH| • In the latter case we 
use the radius of the orbit r' = 1.5 X 10 15 cm = 1/t^o and the magnetic field 
to obtain the Lorentz factor. Finally we consider the GRB021206 (z « 1.25). 
According to Ref. this object has a very compact core with a radius of the 
order of 1 km and a magnetic field « 10 12 Gauss. The synchrotron emission 
region is about 10 8 km from the core so that we estimate the magnetic 
field to be 10 4 Gauss using the transport law Br = const.. From Ref. 92 
we take the cut-off frequency to be uj c = 1 MeV. 

As indicated in Table , the radiation of interest is dominated by very 
high harmonics 10 15 < m < 10 30 exhibiting also large ratios of m/7, typically 
in the range 10 10 < m/7 < 10 22 . The corresponding values for 7 imply also 
(3 ps 1. In this way, the high harmonics present in the synchrotron radiation 
of these astrophysical sources together with the values of the 7 factors of the 
radiating charges highlight the relevance of the large m and large 7 limit, with 
the constraint (m/7) 2 3> 1, to study the induced Lorentz violating effects. 



In a similar way to the standard case j83j . we obtain 

g 2 m^ 1 [ f°° / 3 \ 2/3 

PXm = 777 n T~ JT\ V W dX ^ K ^/3 ( / 

V37T.R 1 + n 2 (Xz m ) I Jm/rhxa V 2m A, ' 
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2m Ac 



4/3 



K V3 I ™ 



for the integrated power in the m th harmonic. Here K p / q denote the Mac- 
donald functions (Bessel functions of fractional order). A first consequence of 
this approximation is the appearance of the cutoff frequency lu\ c — rh\ c loq 
with 

m Xc = ^[l~l3 2 (E)nl(z m )]- 3/2 . (135) 

This name arises because for m > rh\ c the total power decreases as 

P\ m ~ exp(-m/m Ac ). (136) 

Within the same large-m approximation, the integrated power in the m-th 
harmonic can be expanded to second order in £ yielding 



Pm = 
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q z LU 



\/37Ti?7 2 



m c ( m \ 2 „ ( m 
m \m c J Y \m c 



t; muj(3 
7 . 



K 



2/3 



(137) 



where m c = 37 3 /2 and n(x) = x dy K 5 / 3 (y) is the so called bremsstrah- 
lung function. 

Let us notice the appearance of the combination £ ujm /"/ = Mp ) (m /7) 
as the expansion parameter in l|137|l . Here we take M = Mp. As can be seen 
from Table |21 this is not necessarily a small number, which signals the possi- 
bility that such corrections might be relevant in setting bounds upon £. This 
rather unexpected effect is due to the amplifying factor 771/7. Similar results 
have been obtained in calculations of the synchrotron radiation spectra in 
the context of non-commutative electrodynamics |93| . 

Another possibility for observable effects due to £ is to look at the averaged 
degree of circular polarization 



(p + m + p_H)' 



(138) 



where P\ (lo) is the total power distribution per unit frequency and polariza- 
tion A, so that P\(to) = Pmx/ojQ. The average here is calculated with respect 
to an energy distribution of the relativistic electrons, which we take to be 
N(E)dE = CE~ p dE, in some energy range E\ < E < E 2 , chosen as E\ = 
and E2 — > 00 for simplicity. The result is 



(139) 
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( P -3)(3^i)( P + i)r(g + i)r(g + A) 
(P) 3(3p-7) (p-i)r(f + g)^(f + g)' P> ' M (i4Uj 

Again, we have the presence of an amplifying factor in Eq. (|139|1 . given by 
(nuj/qB), which is independent of the form of 77 (p) and not necessarily a 
small number. An estimation of this factor in the zeroth-order approximation 
(£ = = rj), which is appropriate in Ijl39|l . yields (fiu/qB) = u>/ {ujqj) =771/7. 
The expression (I139f) is analogous to the well-known average of the degree 
of linear polarization IIljn — (p + l)/(p + 7/3), under the same energy 
distribution for the electrons [9*4] . 

Finally in Table 0] we have estimated the contributions to the different 
pieces in the phase (I116|l of the Green function, for each astrophysical object. 
Here we have set £ = 1 and M = Mp. The extreme case is 



Table 4. The far field approximation 



Object 


r jr 


uj c /Mp 


(r'/r)((j c /Mp) 


(r'/r) 2 


CRAB 


10" 6 


w -20 


10" 


26 


lO" 12 


Mkn 501 


10- 11 


lO" 15 


10' 


-26 


1Q -22 


GRB 021006 


1Q -24 


1Q -22 


10" 


-46 


1Q -48 



r 

— < 

r 



fa 



< 



(141) 



where all the dependence on £ is negligible in the phase, which reduces to 

(Xz)cj\r-r'\~uj(r-h-r'). (142) 

This corresponds to the CRAB nebulae case, where the assumptions made in 
Ref. |15) are readily verified. A detailed discussion of synchrotron radiation 
in the Myers-Pospelov model can be found in the third reference |52) . 



4.4 General point of view of LIV radiation 

Three paradigmatic examples of Lorentz violating electrodynamics are given 
by the effective theories proposed by Gambini and Pullin (GP) [HOj, Ellis et 
al. (EMN) [2H1, and Myers and Pospelov (MP) [H|. They can be written in 
the general form of Maxwell equations 

V ■ D = 47T/9, V ■ B = 0, (143) 

dB dT> , . 

V x E = - — , V x H - — + 47rj, (144) 

with corresponding constitutive relations 
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D = D(E, B), H = H(E, B), (145) 

which we next present in detail for each case, after reviewing the correspond- 
ing equations. Let us recall that the above equations l|143fl and l|144|) imply 
charge conservation dp/dt +V-j = 0, independently of the constitutive equa- 
tions Ijl45|l • In an abuse of notation we have denoted by £ the electromagnetic 
LIV parameter for all models in the sequel. 

Gambini-Pullin Electrodynamics 

The Maxwell equations for this case are 

V-B = 0, V x (S + 2£V x £ ) + ^ = 0, (146) 

V-£ = 4tt ( o, V X (fi + 2fV X - — =47rj, (147) 

where the electric and magnetic fields are identified from the homogeneous 
equation as 

E = £ + 2£V x £, B = B. (148) 
From the inhomogeneous equations we obtain 

D = £ , H = B + 2|V X B, (149) 

which together with the constitutive relations 

D + 2|V x D = E, H = B + 2£V x B, (150) 

leave the equations in the required form. These equations define the corre- 
sponding functions stated in i|145[l . In momentum space we have 

D = — f E - 2ifk x E + 4£ 2 (k • E) k), H = B + 2i£k x B. 

1 + 4£ 2 fc 2 V / 

(151) 

The admixture of vectors and axial vectors in the constitutive relations pre- 
cludes the parity violation exhibited by the model, together with the presence 
of birefringence. 

Ellis et al. Electrodynamics 

In this case the modified Maxwell equations are 

r)B 

V-B = 0, VxE+ — = 0, (152) 

at 

<9E 

V-E + u- — =^ P eff = 47r(p-u-j), (153) 

V x B - (1 - u 2 ) — + u x — + (u • V) E = 4^j eff 

= 47r(j + u(p-u-j)), (154) 
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which can be written in the form 1(1431144(1 via the constitutive relations 28 

H = B /(w)k x E, D = (1 - f 2 (uj)k 2 ) E + / V)k (k • E) - f(w)k x B, 

(155) 

where we have assumed that u = /(w)k in momentum space. Taking u as a 
vector, this model conserves parity and shows no birefringence. 

Myers-Pospelov Electrodynamics 

This case corresponds to the equations 1(1020 ■ From the last one we can infer 
the constitutive relations 

H = B |d E, D = E-£d B, (156) 

which produce 

V-E=V-D, (157) 

leaving the third Eq. 1(102(1 in desired form. Similarly to the GP case, this 
model violates parity. In momentum space Eqs. ((156(1 become 

H = B + i£wE, D = E + i£wB. (158) 

The above constitutive relations in the three representative models involve 
linear relations among the fields and can be summarized, in momentum space, 
as the local relations 

A(w,k) = ai 3 -(o;,k)£i(a;,k) +py(w,k)-B,'(w,k), 

ffi(w.k) =/3 tf (w,k).B i (w,k) +a ij (uj,-k)E j {u,k), (159) 

where the corresponding momentum dependent coefficients can be read from 
the equations l(151|) . ((155(1 . and ((158(1 . Equations 1(159(1 are the most general 
expressions in which any pair of linear constitutive relations can be ultimately 
written, which allow to express the fields D, H in terms of E, B. 

Parameterization of the constitutive relations 

Let us consider Maxwell equations in momentum space 

k • B (w, k) = 0, kxE(w,k) = wB (w,k), (160) 

ik • D (u, k) = 47rp (u>, k) , ikxH (uj, k) = -iwD (uj, k) + 47rj (u>, k) . 

(161) 

Here we discuss the vacuum situation where the non trivial constitutive re- 
lations arise because of LIV effects. Let us take into account corrections up 
to second order in the LIV parameter £ and let us assume that we are in 
a Lorentz frame where we demand invariance under rotations. This would 
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correspond to the rest frame = (1,0) in the MP model, for example. We 
can always go to an arbitrary frame by means of a passive (observer) Lorentz 
transformation. In this way we have the general expressions 

aiij = a 5i.j + iai&irjkr + hkj, p t j = + ipi&irjh + P2I hkj, 



aoSij + iai£ei r jkr + ct 2 £ hkj, 

(162) 



where a a , Pa , PA , a A , A = 0,1,2, are scalar functions depending only upon 
u>, k — |k|, and £. The property k • B = sets (3 2 — P2 = effectively. In 
vector notation we then have 

D = (a + a 2 k 2 i 2 ) E + ( p a + ic^wf) B + i(p 1 - ia 2 u£) £ k x B, 



H 



o-o + o- 2 k 2 £ 2 ) E + ( Pq + ia^i )B+i(p 1 - ia 2 oj£ l^kxB 



(ft 



(163) 



where we have used the second Eq. (|160fl together with (k'E)k = w(kxB) + 
fc 2 E. 

The generalized Maxwell equations 

Next we substitute (I163|l in Eqs. (|161fl to obtain the corresponding equations 
for E and B. The result is 



i (a + a 2 k 2 ^ (k • E) = 4?rp, 

( 



(164) 



i [a + a 2 k r£ uj~E + i 



k x B 



A + i (^l + Pi) w £ + oi 2 i 2 uJ 2 
+i ((T + Pa) u + i (oliuj 2 - (3ik 2 ) £ B = 4?rj (w, k) . 

(165) 

It is convenient to rewrite the inhomogeneous equations in the compact form 



iP(k ■ E) = 47T/9, 
iojPE + iQk x B + RB = 4irj (w,k) . 



(166) 
(167) 



by defining 



P = a + «2C 2 fc 2 , Q = A + i (f i + Pi) w £ + a2?w 2 , 

i?= (/3 1 fc 2 -a 1 cj 2 )^+i((io + po)^- (168) 

Now we have only three independent functions which depend on to and k. 
The inhomogeneous equations (|166|l and l|167fl can be solved by introducing 
the standard electromagnetic potentials in the Coulomb gauge. The radiation 
approximation is described in terms of the vector potential A x only, which 



Dynamical corrections from space granularity 43 



will exhibit polarizations A = ±1 in the general case and that satisfies the 
equation 

[Qk 2 - Pu? + Xk R] A A = 47rj£. (169) 

Here j^. is the transverse part of the polarized current. The Green function 
G x (uj,k) corresponding to the operator in the LHS of I|ltj9fl satisfies in gen- 
eral a cubic equation in k, which determines the dispersion relations for the 
propagating photon. This means that we have three poles k\(tu). In fact one 
of the poles is due to the (3\ factor in the expression for R in Eq. (|168|l 
and it is located at k ~ Therefore its contribution can be neglected in 
our effective theory valid for k <C f . The remaining two poles will deviate 
little from their position when £ = and will allow us to characterize the 
propagation mode corresponding to each polarization A by a refraction index 
n\{ui) to be read from the appropriate dispersion relation in such a way that 
n\(ui) = k\(u>)/u). The general form of the polarized Green function will be 

G x (cu,r) = ^—F x {uj)e iujn ^ uj)r . (170) 

Notice that from the birefringent case we can go to the non-birefringent one 
by taking n+(uj) = n_(cj) = n(uj), in which case F + (u>) = F~(lu) = F(u>). 

In the following we make explicit some general properties that must be 
satisfied by any electrodynamics characterized by a Green function of the type 
<|170|) . The reality of the electric and magnetic fields leads to the condition 

[G+(w,r)]* = G-(-w,r). (171) 

This implies the relations 

n;H=n_(-w), [F+i^Y^F-i-uj). (172) 

For a birefringent medium the real and imaginary parts of the refraction 
index for circular polarization components can contain both w-even and u>- 
odd terms, provided that they satisfy the first Eq. (|172|l . In the case of a 
non-birefringent medium the real part of the refraction index must be even 
in u>, while the imaginary part must be odd. We can see that the refraction 
indices for the Myers-Pospelov theory, Eq. I)115[l . satisfy these requirements. 
A detailed discussion of this generalized point of view is under preparation 
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